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Abstract 

We present a new mechanism for generation of large-scale magnetic field by thermal convection 
which does not involve the a-effect. We consider weakly nonlinear perturbations of space-periodic 
steady convective magnetic dynamos in a rotating layer that were identified in our previous work. 
The perturbations have a spatial scale in the horizontal direction that is much larger than the 
period of the perturbed convective magnetohydrodynamic state. Following the formalism of the 
multiscale stability theory, we have derived the system of amplitude equations governing the evo¬ 
lution of the leading terms in the expansion of the perturbations in power series in the scale ratio. 
This asymptotic analysis is more involved than in the cases considered earlier, because the kernel 
of the operator of linearisation has zero-mean neutral modes whose origin lies in the spatial in¬ 
variance of the perturbed regime, the operator reduced on the generalised kernel has two Jordan 
normal form blocks of size two, and simplifying symmetries of the perturbed state are now missing. 
Numerical results for the amplitude equations show that a large-scale perturbation, periodic in 
slow horizontal variable, either converges to a short-scale neutral stability mode with amplitudes 
tending to constant values, or it blows up at a finite slow time. 


1. Introduction 


A fundamental problem in astrophysics is to understand the sources of magnetic helds that are 
featured by many astrophysical bodies such as the Sun and the Earth. It is generally accepted that 
the magnetic fields are generated by hydromagnetic processes in the melted or fluid-like interiors 
of the bodies. This idea goes back to J. Larmor [l^ 16|. It is widely believed that the so-called 
o-effect plays a prominent role in such processes. Such a mechanism of magnetic held generation 


was first suggested by E. Parker [l^. It relies on the frozenness of magnetic held into the 


conducting medium, when magnetic diffusion is negligible, implying that a small eddy in turbulent 
flow deforms the magnetic held force line into a loop. If this effect does not disappear when averaged 
over many loops, it gives rise to a mean electromotive force (e.m.f.), which can be parallel to the 
mean unperturbed magnetic held and can amplify the original mean held. 

The theory of mean-held electrodynamics (MFE) is developed around a similar central idea 


2a, ll3|: Suppose the flow and magnetic held are split into the mean and fluctuating parts, and the 


interaction of the fluctuating parts of the flow and the held yields a non-zero mean e.m.f. The MFE 
theory postulates that, at least in the case of magnetohydrodynamic (MHD) turbulence, the latter 
is related to the mean magnetic held (the coefficient of proportionality is traditionally denoted by 
a, and, accordingly, the phenomenon is called the “a-effect”) combined, in some MFE models, with 
the spatial gradient of the mean field. 

The MFE theory does not fully justify such relations. An insight into their mathematical roots 
is provided by the multiscale stability theory (MST). It treats an idealised case, in which turbu- 














lence is modelled by a small-scale laminar flow, periodic in space and time (see 33|]). In contrast 


with MFE, mathematically rigorous results are then obtained by applying the asymptotic theory 
of PDE homogenisation to the problem of linear or weakly nonlinear stability of small-scale states. 
Purely hydrodynamic UM perturbations, the kinematic dynamo problem B [ 35 I. Bl (which is 
an instance of the general linear MHD stability problem, where the perturbed state is amagnetic, 
the flow and magnetic components of the perturbation therefore decouple, and one focuses on the 
magnetic perturbation), and full perturbations of forced MHD or convective MHD states [Il,& 
3l| (see also j^. Chap. 6-9) were considered. Perturbations are supposed to involve spatial and 


temporal scales that are much larger than the respective periods of the perturbed states. The per¬ 
turbations are linear combinations of amplitude-modulated small-scale modes (i.e., eigenfunctions 
of the linearisation whose periods coincide with those of the perturbed states) associated with the 
same eigenvalue. The coefficients depend exclusively on slow temporal and spatial variables. They 
are usually called amplitudes, and the equations governing their dynamics are called amplitude 
equations. The evolution of such large-scale perturbations is essentially controlled by the asso¬ 
ciated eigenvalue of the constituting small-scale modes — only neutral (belonging to the kernel 
of linearisation) small-scale modes can instigate instability in the presence of large scales; conse¬ 
quently, MST usually focuses on large-scale amplitude modulation of neutral modes. Since typically 
small-scale neutral modes have non-zero means, amplitude equations for large-scale perturbations 
of forced MHD states (convective or not) are mean-held equations similar to those considered in 
MFE theory. However, as we will see, this is not always the case: translation-invariant physical 
systems, such as free thermal hydromagnetic convection in a horizontal layer investigated in the 
present paper, possess zero-mean neutral small-scale stability modes. Consequently, the mean-field 
description of the dynamics of large-scale perturbations of such physical systems is inadequate. 

We intend to carry out a detailed investigation of stability to large-scale weakly nonlinear per¬ 
turbations of the steady and time-periodic regimes of magnetic held generation by free thermal 
convection of rotating electrically conducting fluid in a horizontal layer, that were determined 
numerically in and (iil. Here we consider space-periodic convective MHD steady states con¬ 
stituting the branch [7|. The group of symmetries of this steady state is generated by the 
symmetry about the vertical axis, X 3 , (note that this is not axisymmetry, see section [2]) [3, 33| and 
the superposition of reflection about the midplane with translation by half a period in the hori¬ 
zontal direction xi. This group of symmetries is smaller than other ones, typical for the parameter 
values considered ibid., and the system of amplitude equations, derived in j33(|, is inapplicable for 
large-scale perturbations of states having this group. The symmetry about a vertical axis implies 
that the steady state does not possess the a-effect. The large-scale dynamics is due to the interplay 
of the combined eddy diffusivity and eddy advection. We And that their interaction cannot sustain 
stationary generation of large-scale magnetic held: a large-scale perturbation from the class which 
is described by the multiscale formalism either converges to a short-scale neutral stability mode, 
or it blows up at a finite time. 

The paper is organised as follows. In section [2] we derive the system of amplitude equations for 
steady states of free hydromagnetic convection, which have the same group of symmetries as those 
comprising the branch In section [3] we present results of numerical analysis of the system of 
amplitude equations for several states belonging to this branch. Finally, we make remarks triggered 
by our investigation. 


2. The multiscale formalism for large-scale perturbations 
of free hydromagnetic convection 

In this section we derive amplitude equations for large-scale perturbations of small-scale steady 
free hydromagnetic convection in a horizontal layer of electrically conducting fluid rotating about 
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the vertical axis. A field, depending only on the fast spatial variables, x, and time, t, is called 
small-scale; if, in addition, the field depends on the slow horizontal spatial variables, X = e{xi,X 2 ) 
and on the slow time, T = e^t, where s > 1, it is called large-scale. We will use asymptotic methods 
that are standard within the MST approach. We assume that the perturbed state is symmetric 
about the vertical axis and has equal periods ^ in xi and X 2 ; further assumptions are introduced 
where they become relevant. 

A three-dimensional field f is called symmetric about a vertical axis passing through point 
(ai,a 2 , 0 ) when the following conditions are satisfied 0 , 0 : 


/i(ai -xi,a2 -X2,X3) = -/i(ai + Xi,a2 + X2,X3), 
/2(ai - Xi, 02 - X2, X3) = -/2(ai + Xi,a2 + X2, X3), 
fsiai - Xi, 02 - X2, X3) = /3(Oi -h Xi, O2 + X2, X3). 


In particular, this symmetry implies that the flow is vertical everywhere on the axis. Together with 
the ,^-periodicity in the horizontal directions, the symmetry about the Cartesian axis X 3 implies 
that the field is also symmetric about vertical axes through points <^n/2, where n = (n,i,n 2 , 0 ) 
has integer components. 


2.1. Small-scale convective hydromagnetic steady states and their perturbations 

The state, whose stability we examine, is governed by the Navier-Stokes, magnetic induction 
and heat transfer equations. In the coordinate system, co-rotating with the fluid layer about the 
axis X 3 , they are 

dY/dt = z/V^V + V X (V X V) - H X (V X H) + /^©eg + rV x eg - VP, (1) 

d'B./dt = + V x (V x H), (2) 

dQ/dt = kV^© - (V ■ V)© + 51 / 3 , (3) 

V ■ V = V • H = 0. (4) 

Here V is the flow velocity, H magnetic field, Q = d — di + 6 x 3 the difference between the 
temperature, il, and the steady state linear profile, ili and '^2 prescribed temperatures at the upper 
and lower horizontal boundaries xg = 0 and 1, respectively, 5 = "di — '^2 > 0 and r = \/Th. The 
Navier-Stokes equation ([1]) involves buoyancy, Coriolis and Lorentz forces. No external body forces, 
heat or electric sources are present, i.e., there are no source terms in equations ([I])-(|l]); in this case, 
we call convection free. Since the small-scale state W = (V, H, ©) is supposed to be steady, the 
time derivatives in the l.h.s. of ([I])-® vanish. 

Equations ([I])-(I1]) are invariant with respect to the position of the axis of rotation, i.e., the 
equations are the same for any vertical axis of rotation. In an inertial (e.g., steady) coordinate 
system, a periodic arrangement of rotating convective cells, with the pressure periodic in horizontal 
directions, is impossible. Periodicity in xi and X 2 becomes possible after the pressure is corrected by 
r^p^/ 8 , where p is the distance to the axis of rotation (see j^); this term represents the centripetal 
force preserving the arrangement of cells (regarded as rigid bodies) rotating with constant angular 
velocity r/2. Clearly, rotation about a vertical axis is compatible with the symmetry about a 
vertical line; however, since the governing equations are invariant with respect to the position of 
the axis of rotation, the axis of rotation does not necessarily pass through the centre of a periodicity 
cell or coincides with the axis about which the flow is symmetric. 

We consider horizontal boundaries of the fluid layer that are stress-free, ideally electrically 
conducting and kept at constant temperatures. This translates into the following equations: 

dVi / dx ^ = 91/2/9x3 = 1/3 = 0, 

dHi/dxs = dH2/dx3 = ifg = 0, 

© = 0 . 
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(5) 

( 6 ) 
(7) 



We define the mean over the fast spatial variables, (•), and the flnctnating, {•}, parts of a held /: 

(/(X,T,x,f))^— / / / /(X,T,x,f)dx, {/}^/-(/). 

Same dehnitions hold for a vector held f, and we also dehne the mean of its horizontal part, 
(f)/j = (/i)ei + (/ 2 )e 2 , and the hnctnating part {f}^ = f — (f)/j. Later (see section [2.2.ip it will 
become clear why these averaging procednres agree with the mathematical natnre of the problem. 
We note that the Coriolis force can be reexpressed as t{V}^ x 63, since the constant horizontal 
force T{V)h X 63 is a gradient of a linear fnnction of horizontal coordinates that modihes pressure. 

Upon this modihcation, the operator of linearisation of equations ([I])-([2D around the steady 
state W, which we denote by il = (£*', 2,^), takes the form 

£’'(v, h, 0) = - dv/dt + + V X (V X v) + V X (V X V) 

— H X (V X h) — h X (V X H) + /^des + r {v}^ x 63 — Vp, 

£'*(v, h, 0) = - dh/dt + r/V^h + Vx(Vxh + vxH), 

£®(v, h, 0) = - dO/dt + - (V ■ V)e - (v ■ V)0 + Sv 3 . 

We consider a perturbation of the steady state W, whose amplitude is 0(e). The perturbed 


state, (V + ev, H + eh, 0 + e9), satishes equations ([I])-(|l]), whereby 

i:''(v, h, 0) = - e(v X (V X v) - h X (V X h)) + Vp, (8) 

i:^(v,h,0) =-eV X (v X h), (9) 

£®(v,h, 0 ) =e(v-V) 0 , ( 10 ) 

V ■ V = V ■ h = 0. (11) 


Following the recipes of the homogenisation techniques, we now assume that the perturbations 
depend on the slow horizontal spatial variables X = e(xi, 0 : 2 ) and on the slow time, T = e®t, where 
s > 1 , and expand the perturbation in a power series in the small scale ratio e: 

00 

V = ^v„(X,T, x,t)£”, 

n=0 
00 

h = ^h„(X,T,x,^)£^ 

n=0 
00 

0=5^0„(X,T,x,t)£T 

n=0 

Upon substituting these series into the governing equations B-imii, we obtain a hierarchy of 
systems of equations. In principle, it can be solved at any order, and thus complete series 
can be reconstructed. In this paper we are only interested in the leading-order terms and therefore 
will only consider the hrst 3 systems in the hierarchy. 

The solenoidality conditions flTT]) at order e"' yield equations, whose mean and fluctuating parts 
for any n > 0 are 


Vx ■ {Vn)/i — 0, 

(15) 

Vx ■ {iin)h = 0, 

(16) 

Vx • v„ Vx ■ = 0, 

(17) 

Vx ■ h„ -|- Vx ■ = 0 

(18) 


( 12 ) 

(13) 

(14) 
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(any term of expansions with a negative index is zero by definition). We can now comment on the 
modihed form of the Coriolis force term in the operator in (|H]): By virtne of solenoidality 
of the mean held, we can introdnce a stream fnnction, ijj, whereby (v)/i = {—d'ip/dX 2 , d'll^/dXi, 0 ). 
Thns, the difference between the original, rv x ea, and modihed, t{v}^ x 63 , forms is eqnal to 
r(v)/i X 63 = rVx'05 which is compensated by modifying the pressnre pertnrbation, p. So, althongh 
now the how pertnrbation, v, depends on the slow spatial variables, the operator of linearisation, 
of the Navier-Stokes eqnation ([T]) can still be nsed in the form dH]). 

We assnme henceforth that the pressnre pertnrbation p is i^-periodic in the fast horizontal 
variables. 

2.2. Orders^ equations 

This section is devoted to the stndy of the generalised kernel of the operator of linearisation £. 
We hnd that the operator acting on the generalised kernel has two Jordan normal form blocks of 
size two, and, conseqnently, the means of the how and magnetic components of nentral modes of £, 
Sfc for k = 1,2, are related by a matrix U of size two introdnced below. The pertnrbed states are 
missing the symmetries which make this matrix vanish. As a resnlt, it will tnrn ont that becanse of 
the solenoidality in the slow variable of the spatially averaged pertnrbations of the form ffT^ - fflT)) . 
the latter can be described by a self-consistent system of amplitnde eqnations only if they depend 
on a single spatial variable Y = q • X, where a nnit vector q is normal to an eigenvector of U. 

2.2.1. The Jordan normal form of the operator of linearisation 

in the invariant subspace associated with the zero eigenvalue 
The hrst system of eqnations in the hierarchy is obtained at order 

i:(vo, bo, 6»o) = 0, Vx ■ Vo = Vx ■ bo = 0, (19) 

i.e., the held (vo, bo, 6 * 0 ) belongs to the kernel of the operator of linearisation, £; it mnst satisfy the 
reqnirements for pertnrbation, i.e., the solenoidality conditions for v and h, .^-periodicity in the 
fast horizontal variables and bonndary conditions ©-©■ 

As nsnal, we derive the adjoint operator £* = (£*’', £*®) in the Lebesgne space L 2 of 

solenoidal vector helds by nsing vector analysis identities and, essentially, integrating by parts the 
l.h.s. of the dehning eqnation 

i 2 w ■ w* = w ■ : ( 20 ) 

£*^(v, h, 9) = dv/dt + z/V^v - V x (V x v) 

+ ^(H X (V X h) - V X (V X V) + sees - r {v} x 63 + 0V0), 

£*^(v, h, 9) = dh/dt + pX^h + V x (H x v) + q3(v x (V x H) - V x (V x h)), 

£*®(v, h, 9) = 89/dt + kX‘^9 + (V ■ X)9 + fiv^. 

Here ^ is the projection onto the snbspace of solenoidal vector helds satisfying the reqnirements 
for pertnrbation. 

We assnme henceforth that the small-scale pertnrbed steady state W = (V, H, 0) is stable to 
small-scale pertnrbations, i.e., real parts of all eigenvalnes of £ are non-positive. It is then natnral 
to assnme that the evolntion of large-scale pertnrbations of W is examined after all the fast-time 
transients have decayed. Since the pertnrbed state is steady, we will assnme that the pertnrbations 
do not depend on the fast time, and hence the derivatives in t in the dehnitions of the operators 
£ and £* will be omitted. (It wonld be necessary to preserve them when investigating large-scale 
stability of time-periodic states.) 


5 


Since £*(efc,0,0) = £*(0,efc,0) = 0 for fc = 1,2, the kernel of the operator £ is at least four- 
dimensional. Differentiating ([I])-(ll]) in Xk for k = 1,2, we hnd that 

Sfc+2 = dW/dxk (21) 

belong to the kernel of £. Furthermore, 

£(efc, 0, 0) = -Sfc +2 e ker £; (22) 

therefore, the Jordan normal form of £ involves two blocks of size 2 associated with the eigenvalue 
zero, and £ has at least 6 eigenhelds associated with this eigenvalue, at least two of which are 
generalised. Henceforth, we assume the generic case where zero is a 6-fold (generalised) eigenvalue 
of £. Clearly, {£,^{w))h = (£^(w))/i = 0 for any w = (v, h, 0) satisfying the requirements for 
perturbation. Consequently, it is simple to show that any eigenheld, a generalised one or not, which 
has a non-zero mean horizontal part of the flow and/or magnetic component, is associated with the 
zero eigenvalue. Such eigenhelds exist, because eigenhelds of the elliptic operator £ constitute a 
complete basis in the Lebesgue space L 2 ; hence, the kernel of £ involves the helds = (S^, S^, S^) 
for fc = 1, 2, satisfying the requirements for perturbation and such that 

2 

£St = 0, = {sa^ + ^ sj = {sq, + e,. 

m=l 


The first auxiliary problem is to hnd the helds and the matrix U = 


Since 


Mil ^12 

■^21 U 22 

normal Jordan forms of an operator and its adjoint coincide, there exist two generalised eigenhelds 
S* = (S™, S*^, 5*^), n = 1 , 2 , associated with the eigenvalue zero of £*, i.e., such that 


£-s:= 0 • (s:’a=(sa)»= 0 , 


v,.sa=v,.s:'‘=o 


(23) 


y. m=l 


and also satisfying the requirements for perturbation. It is straightforward to demonstrate that 

wji u '^2 

U 21 U 22 

where the superscript t denotes the transpose of a matrix. 

By linearity of the system of order equations flT^ . its solution is 


= -{U 


t\-i 


(24) 


Wo = (vo, ho, ^o) = ^ C 0 fc(X, T) Sfc(x), 
k=l 


Po = 5^c„t(x,r)sr(x), 

k=l 


(25) 


where po and denote the gradient potentials arising in the equations £^wo = 0 and £^Sfc = 0 , 
respectively. The horizontal mean of the magnetic component of the hrst relation in fl25|l reduces 
to cofc = {ho,k) for k = 1,2, and hence {vo)h = U{ho)h- 

The solvability condition for a problem £w = (P, f^, /^) consists of orthogonality of the r.h.s. to 
the kernel of the adjoint £*. Since the kernel of £* is comprised of constant vectors in the horizontal 
flow and magnetic components, the solvability condition reduces to (P)^ = {f^)h = 0 . It is 
straightforward to establish that this condition is necessary; that it is sufficient follows from the 
Fredholm alternative theorem applied to the problem (V^)“^£w = (V^)“^(P, f^,/®) (the proof is 
essentially the same as in section 3.2 in j^). It is due to this form of the solvability condition 
that spatial averaging over the periodicity domain is the only appropriate kind of averaging in the 
problem under investigation. 
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2.2.2. Solenoidality of the mean fields 

Here we discuss how to satisfy the solenoidality conditions flT^ - flThl) . We assume that ampli¬ 
tudes cofc are space-periodic in the slow variables X = (Xi, X 2 ) and expand the amplitudes in 
Fourier series, which it is convenient to express in the form 

(vo)fe = ^ (vpfcqei-h n 2 ,fcqe 2 ) (ho)h = ^ + h 2 ,kqe 2 ) (26) 

fe, |q| = l fc,|q| = l 


solenoidality of 


is equivalent to the orthogonality condition Wq-*- ■ q = 0. 


and consider the conditions independently for each unit wave vector q = (gi, 52 )- 

Solenoidality of (ho)h then requires (hi ^q, h 2 ,A:q) = C'^qH'*', where q-*- = (^ 2 ,—<?i), and hence 

^Ijfcq _ ^ h\kc{ 

^2,fcq h2jfcq 

It implies q-*- || Ucfi-] thus, q-*- is an eigenvector of Li. We denote the associated eigenvalue by 
A = q-*- -Ucfi-. In the basis {q‘'',q}, the matrix U is upper triangular; consequently, the wave 
vector q is an eigenvector of lA^ associated with the second eigenvalue A' = q ■ Wq. 

Thus, by virtue of the solenoidality conditions the series fl26|l may involve only two wave vectors. 


q^ and q^, that are normal to eigenvectors of W; hence fvo)h and (ho)h split into sums of vector 
helds, each depending on just one spatial variable, W = qW X or F 2 = ■ X. We will see 


in section [231 that the system of amplitude equations for cqa, is nonlinear (cubic). Therefore, a 
weakly nonlinear perturbation depends on one spatial variable H = q ■ X such that q-*- is an 
eigenvector of U (otherwise, products of functions of Yi and 12 are generated by the nonlinearity 
of the governing equations for Cok)- Perturbations with initial amplitudes that violate this one- 
dimensionality restriction cannot thus be expanded in the series flT 2 |) - flTT)) and are not described 
by homogenised equations that we will derive. Arising for L( ^ 0, this is the main difference with 
the cases considered in Summarising, 


(^^0,1) 

(^^0,2) 


= XC{Y,T)q^, 


(ho,i) 

{ho,2) 


= C{Y, T)q^, Wq^ = Aq^, F = q ■ X. 


(27) 


We will henceforth use the same notation for the two-dimensional vector q and the three- 
dimensional horizontal vector, whose horizontal component coincides with q; similarly for q-*-. 


2.3. Orders^ equations 

At order we obtain the system of equations 


'2*^(vi, hi, ^i) -|- 2i/(Vx ■ Vx)vo -|- V X (Vx X vq) — H x (Vx x ho) 

+ Vo X (Vx X Vo) - ho X (Vx X ho) - VxPo = 0, (28) 

i:^(vi, hi, ^i) -F 2?7 (Vx ■ Vx)ho -1 Vx x (vo x H -F V x ho) -1 Vx x (vo x ho) = 0, (29) 

i2^(vi, hi, efi + 2«(Vx ■ Vx)0o - (V ■ VxWo - (vo ■ Vx)0o = 0. (30) 


For the slow time T = et, these equations would also involve the derivatives of the respective fields 
in the slow time. The solvability condition then defines the operator of the significant combined 
a-effect: it consists of spatial means of linear (in wq) terms in the horizontal components of the 
flow and magnetic equations. However, the perturbed state W = (V, H, 0) is supposed to be 
symmetric about the vertical axis, X 3 . Consequently, solutions to the hrst auxiliary problems 
are antisymmetric, the r.h.s. of the equations are symmetric, and the a-effect operator is zero: 
convective hydromagnetic states with such a symmetry lack any signihcant (i.e., responsible for 
the evolution of large-scale amplitudes) a-effect. In the absence of signihcant a-effect, we have to 
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switch to the slow time T = e'^t. Since the solvability condition is verihed, a solution to this system 
exists and, by linearity, takes the form 


Vl 


hi I ^ ^ j CifcSfc ^ ^ Gmk 


dct 


Ok 


where G 


mk I 


V ) 

' k=l \ 

t) for k = 1, 

..., 4 and m - 

_ 

^ ^mk — 

2 vdS1/dxm 

= - 

2 r]dSk/dxm 

_ 

^TTlk 

2 KdSl/dxm 


m=l 


dX„ 


^ ^ Qm/cCo/cCom j 5 


771=1 




V ■ or,.* = - , V . Gl. = - {S*„} 

and Qmfc(x, t) for m,k = 1, ...,4 are solutions to the third auxiliary problem: 


(31) 


(32) 

(33) 

(34) 

(35) 


2£’'Q„i = - SJ X (V X S:;.) + S'; X (V X S") - X (V X SI) + SO. X (V X SO), (36) 

2£'‘Q„i = -Vx(SOxS0. + SO,xS0), (37) 

2£'’Q„i = (SO ■ V)S0. + (SO.. V)Sl (38) 

V-Q0.i=V-Q0.i = O (39) 


(for the symmetrised last sum in fl3Tl) h The potential pi in fl28|) can be expressed in terms of the 
gradient potentials and arising in fl52D and fl5BD . where the flow components and 
satisfy relations fl5SD and fl5^ . respectively: 


pi = Yl 


GkSl + ^ 


711=1 


dcpk 

dXm 


E 

771=1 




(Pi) = 0. 


For the same symmetry reasons as for the full equations fl25D - fl5UD . the r.h.s. of the equations in 
formulations of the second and third auxiliary problems are symmetric, and hence the solvability 
conditions for them are satished. Symmetric and antisymmetric helds constitute invariant subspaces 
of the operator of linearisation £, and thus Gmk and Qmk are symmetric about the vertical axis. 

Examining the system of equations for the difference G 14 — G 23 , it is straightforward (albeit 
tedious) to show that 

Gi 4 — G 23 = (V X 63 , H X 63 , 0). (40) 

Differentiating the equation £(Sfc) = 0 in xi and X 2 and comparing the results with fl36|) -fl39 |) , 
we hnd 

dSk/dXm-2 = 2Qmfe = “^Qkm, (41) 

for any /c = 1,..., 4 and m = 3,4. 


2.4- Orders'^ equations 

At order we obtain the system of equations 

dvo/dT = 22 {v2, h2, 62) + 2 z/(Vx ■ Vx)vi + z/VxVq 

+ V X (Vx X Vl) - H X (Vx X hi) + vq x (Vx x vq + Vx x vi) (42) 

- ho X (Vx X ho + Vx X hi) + vi x (Vx x vo) - hi x (Vx x ho) - VxPi, 

dho/dT = £"(v 2 , h 2 , 62 ) + 2i7(Vx ■ Vx)hi + T^V^ho (43) 

+ Vx X (vi X H + V X hi + Vo X ho) + Vx X (vq x hi + vi x ho), 
dOo/dT = £^(v 2, h2, 62 ) + 2k{V^ ■ Vx)0i + 

- (vi ■ Vx)0o - (vo ■ Vx)0i - (vo ■ Vx)0o - (V ■ Vx)0i. 


(44) 




Two equations for the evolution of amplitudes of large-scale perturbations are the solvability 
conditions for the system fl42l) - fl44p : orthogonality to the kernel of the adjoint operator, £*, i.e., 
the requirement of vanishing of the horizontal means of the flow and magnetic equations. By fl42|) . 


^(vo)fe 

dT 



m=l 


d'^cok 

'jmkdX.dX„ 


-Ea; 

m=l 


jrak 


^(^OmCo/c) 

~dX~ 


Vx(Pi)/i, 


(45) 


where we have denoted 


= (V X (e, X G:;,) - H X (e, x 
= (V X X q:;,) - h X (e, X Qt,) - 


and from 


^(ho)fe 

dT 


4 

hVx(ho)/, - eg X VxJ]] 

k=l 


V m=l 


h 

mk 


dc, 


Ok 


dX„ 


m=l 


h ^ 

m/c^Or) 


.^Ok 


(46) 


where 

Ol = ((VxG:;.j-HxG::.J.e3), 

At, = ((V X QD., - H X Q;, + s; X Sj) . e,). 

The quantities and describe eddy diffusivity, and and eddy advection. 

Signihcant simplihcations stem from relations (127)) . The last term in flT5|l serves to eliminate 
from this equation the component, parallel to q; consequently, fl45|l is equivalent to the scalar 
equation 


.dC_ d^C 



D 


jni, 


kQjQr. 


d^CQk j_ 


dY^ 


A 


jm, 


kQj 


^(COmCofc) I 


m=l 


dY 




(47) 


All terms in fl46|) are parallel to q-*-, and, in view of (ITT)) , it reduces to the scalar equation 


ar ^ ’•ar^ ar ' 

k =3 


(48) 


where 


2 2 

a'‘ = EE(-i)’”^‘a;;.,,3-™?3-., 

k=l m=l 


k=l m=l m=l 

Relations fl4D]) imply the following identities for eddy coefficients in (l45]) - flT^ : 

Dt = Dt„ D”,, - D”33 = (-(V X e 3 )n + (H x e^)H,), 

for j = 1,2. Expressions for coefficients in the last sum in (ITT)) can be simplihed by using an 
identity that follows from fITT]) : 


■ (A;„, + = E 

j j 

where m = 3,4 or fc = 3,4. 
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Since the two equations fH7j) and fHS]) involve the slow time derivative of the same quantity, (7, 
a simpler non-evolutionary equation can be derived, equivalent to any of the two. Multiplying (H5]) 
by A, subtracting it from fH7|) and integrating the result in Y we hnd 

Y ( - v>‘)+ ( 49 ) 

\ j=l m=l k=l 

4 2 2 \ 2 4 4 

^OA: | AvA C 

k=3 j=l m=l / j=l k=l m=l 

where K{T) is a function of slow time alone. We will consider perturbations (in particular, cqa,) 
that are periodic in F; by appropriately rescaling e, the period then can be made equal to 27r (in 
agreement with fl26H i. Averaging fH9|) in Y over a period uniquely determines K{T). 

2.4-1- The closing evolutionary equation 

Another evolutionary equation will be derived from fH2]) - fm|) owing to the fact that the gener¬ 
alised kernels of the operators £ and £* are six-dimensional. We begin by making a general remark 
on the means of the perturbation helds, which will be used in the derivations. Let ((•)) denote the 
averaging in Y over a period: 

{f(Y,T)))^T r f{Y,T) AY. 

Clearly, application of the two averagings, (■) and ((•)), to the flow and magnetic equations obtained 
at order yields 9(((v„)))/c}T = 0 (there is no contribution from nonlinear terms by virtue of 
(lT5|l - (lT8ll ) and 9(((h„)))/5T = 0. Thus, the means over both slow and fast variables of each term 
of the expansions of the flow and magnetic perturbation are independent of time. We investigate 
instabilities (if any) of the steady state W developing in their own right and not due to a mean 


flow through the layer, or an imposed mean magnetic held, and therefore demand 

((K))) = (((h„))) = 0. (50) 

We split V 2 and h 2 , into solenoidal and gradient parts: 

V 2 = vf + Vx<, Vx-vf = 0, VM’' = -Vx-Vi, (51) 

h2 = hf + Vxhr, Vx-hf = 0, V2hr = -Vx-hi, (52) 

where we have used relations (ITTD and flTSD for n = 2. Upon substituting these sums into flT2|) -fHi )) , 

we recast this system as 

Qwf = dwo/dT - F, (53) 


where W 2 °* = (v 2 °\ h 2 °\ 62 ) and F denotes the sum of all the remaining terms in these equations. 
Scalar multiplying (l53il by S*, j = 1, 2, we hnd using ([23]): 

^{w„ . S*) - (F . S') = (£wr' • S-) = (wf • £-S*) = ^ uY{v,„} + 

m=l 

(note that, by our constructions in section 12.2.11 relation fl20|) dehning the adjoint operator £* is 
only valid for helds w and w*, whose how and magnetic held components are solenoidal in the fast 
spatial variables — this has forced us to perform the Helmholtz decompositions (l5Tll - (j52ll ). Now, 
by virtue of (El), 

Y (^<"'» ■ s;) - (F ■ s;)) e, = + (h,)i. (54) 

j=i V / 
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However, by flTB]) and ffTBD for j = 2, (v 2 )/i and (h 2 )/i are solenoidal in the slow variable X, and, 
since their spatial dependence is exclusively on H = X-q, they are parallel to q-*- (Fourier expansion 
in Y can be used to show this for zero-mean helds; see fl50l) for n = 2). Furthermore, since q-*- is an 
eigenvector of W, the r.h.s. of fl5^ is parallel to q-*-. Thus, scalar multiplying fl5T|) by q, we obtain 
a closed equation in cok'- 

E(E^(Si.S-)-{F.S-)'),,=0. (65) 

2=1 \k=l J 

It can be simplihed using relations fl22|) and fl23|) : 

(S,+2 ■ Sp = -{^{ek, 0, 0) ■ Sp = -4^.. (56) 

Since q is an eigenvector of the matrix fl2n) . 


^(Sfc+2 ■ S*)qj = - = efc • (W*) ^q = 


By virtue of 


i=i 

takes the form 


i=i 


a 2 4 2 

-giS2) ■S*)C' + —^gfc_2Cofc) = ^^^(F-S*). 


dT\ 




k=3 


(57) 


i=i 


It is straightforward to calculate the scalar products in the r.h.s. of this equation, using the 
expressions for the two leading terms in the expansion of perturbation, that were found in sec¬ 
tions 12.2.11 and 12.31 The result has a simple structure, but involves bulky coefficients: 


2 2 4 


(F • SJ) = E E E + ttt 


72=1 771=1 k=l 


72=1 771=1 k = l 


dX„ 




(58) 


4 4 4 


EEE ^On^Om^Okfi 


Lkr 


72=1 771=1 k = l 


where 

fLtj ^{{‘ivdGlJdXr, + + V X (e„ X G;.,) - H X (e,. x G^j) ■ S“ 

- + {iridGlJdx. + + e„ X (V X G^. - H x G;.,)) ■ Sf 

+ {2^dGlJdxn + - V„Gl,)sf - £(ves.t„ VGU.O) ■ S'), 

flmkj ^{{‘i‘'dQ'UldXn + 2V X (e„ x Q;;^) - 2H x (e„ x Q*,.) 

+ SC. X (V, X Gl, + e„ X sa - Si x (V, x Gj, + e„ x Sj) 

+ Ca X (V, X SCJ - Gl, X (V, X S^)) ■ S™ - 

+ (4i)SQL/&» + e„ X (2V X Ql, - 2H x QC.* + SC. x Sj + SJ x Sj.) 

+ V, X (Sa X Gl, - Si X GCa) • Sf + (4xaQljax„ - 2V„Ql, 

- (GC, ■ V,)Sl - (Sa ■ VJGl - Sl,„stjSf - S.(VQl^„ VQl„.0) .S'), 
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=( (s> (Vx X q:;.j - s; x (Vx x q:;^) + q:, x (Vx x sy 

- QL X (Vx X S^)) . S™ + (Vx X (S” X Ql - S^: X Q;,)) • Sf 

- ((Qr„* ■ Vx)S“ + (S” . v,)Ql,)sf), 

5” is the Kronecker symbol, and Q’^^n are space-periodic solutions to the 

Poisson equations 

_ (~iv 

^ ^mkn ^mk^n'i ^ ^mkn ^mk^n'i ^ ^mkn ^mk^n'i ^ ^mkn ^rak^n' 

The following identities for the coefficients / prove useful: 

1 . For n = 3 or 4, and any m, k and j, straightforward algebra with application of fl2Tl) and fj4Tl) 
yields two relations: 


flmkj + fLmj = 2{{d2.Qmk/dXn-2 “ S.{dQmk / dXn- 2 )) ' S*) 


and 


fmnkj d* fmknj d“ fkmnj d* fknmj -2{{dS.Qmk/dXn-2) ■ Sp. 

Consequently, no terms involving factors Cqs or co 4 are present in the last sum in the r.h.s. of 
2. For any n,j = 1, 2, 


/n34j ~ /) 


nA3j ’ 


(69) 


since by inspection 

f2 _ f2 -/p pG„3 dGnA \ gA_Q 

Derivation of equations for amplitudes is completed by performing the following operations: {i) 
substitute in the r.h.s. of fj4^ and fl58l) the variables cqi and co 2 , using the second relation in fl2711 : 
(ii) in the l.h.s. of fl57|l . eliminate the term, proportional to dC/dT, by subtracting the appropriate 
multiple of flT 8 |l : (in) introduce three new variables: = C, ^2 = QiCos + ? 2 Co 4 (the l.h.s. of the 

modihed flSTl) is proportional to 8 ^ 2 /dT), and .^3 denoting the linear combination of C, Cqs and 
Co 4 , whose derivative in Y is the l.h.s. of fl4^ : {iv) in the modihed fl57D . replace d‘^^s/dY‘^ by a 
derivative of a quadratic form of applying the non-evolutionary equation flT^ . Upon expressing 
the r.h.s. of (HHD . the modihed fISTD and flTOD in terms of we obtain with the use of fl5^ three 
equations. 


dT 

d ^2 

dT 

d^3 

BY 


52 


gY 2 


Y ^ gc2 

+7i^) 


k=l 




gc , 

+ /^22 C 2 ) + ^ 'l2ij 




73ij iiij ^ l3ij 






(60) 

(61) 

(62) 


The coefficients f3, 7 and a can be calculated following transformations {i)-{iv) and using the 
expressions that we have obtained for the coefficients of equations fl45]) . flFT)) and fl4^ in terms of 
solutions to the three auxiliary problems and eigenmodes from the generalised kernel of £*, the 
adjoint to the operator of linearisation. 

We hnd from (l60|) that ((^ 1 )) does not depend on the slow time; hence, in view of (l50|) for n = 0 
and the hrst equation in ([27]), ((^ 1 )) = 0. Relation (162|) does not completely specify ^ 3 , since it 
reduces to a trivial identity for the zero wave-number Fourier harmonic (in Y) of fl62|) . We thus 
lack an equation to close the system (16U]) - (I62D . 
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2.4-2. The closing evolutionary equation 

To obtain the closing relation for the leading-order terms of the pertnrbation, we have to consider 
the eqnations obtained at even higher orders. Snch a sitnation is often enconntered in the dynamo 
theory, for instance, one has to descend by 2 orders of magnitnde to derive the magnetic a-effect 
in the classical almost-axisymmetric Braginsky dynamo j^. 

In order to derive the missing eqnation, we average (l5T|l in F, apply condition (jHil for n = 2, 
and find 

E ■ S*> = (({F» ■ S*). (63) 

k=3 

The component of (j54ll . parallel to q, has been fnlly explored in Section [2.4.1( we consider now the 
component, perpendicnlar to q. Scalar mnltiplication of (l63ll by q-*- and the nse of (j56jl yields 


E (E“t^K*«)) +««■ s]) I = 0. 

j=i \k=i ' 

We express here the mean of 

2 4 4 


(64) 


(((f».s-)=eee((|w<^»”>)) fnmkj T ^ ^ ^ ^ ^ ^ ^^OnCOmCok]) fnmkj y 


2 2 2 


n=l m=l k=l 


n=l m=l k=l 


in terms of nsing fl5^ : 

((hi6 + h'2^3} = ^6^(726 + 73^3)^ + 

and snbtract from this eqnation flbTll averaged in Y and mnltiplied by a' ja. This finally yields the 
closing eqnation: 

o p\ 

((hl6 + /^2^3)) = < 1 ^ 6^(7426 + 743 ^ 3 )^- (65) 

Eqnations fl60|) -fl6 ^ and (l65|) comprise a closed system in ^j, eqnivalent to the system of amplitnde 
eqnations. 


3. Numerical results 


In this section, we present resnlts of nnmerical investigation of the two-scale weakly nonlinear 
stability of space-periodic steady-state convective dynamos in a horizontal layer of rotating con- 
dncting flnid, that comprise the branch labelled in j^. The branch is parameterised by the 
Taylor nnmber, Ta, and exists in the interval 673.7 < Ta < 704. These convective MHD steady 
states are symmetric about the vertical axis, 0 : 3 , and hence the significant (acting on the ampli¬ 
tudes) a-effect is absent. Dynamos from this branch have the group of symmetries D 2 , the second 
group generator being the composition of reflection about the midplane and the shift by half a 
period in the horizontal direction xi. For this group of symmetries, the matrix U is non-zero, this 
invalidating the multiscale analysis of large-scale perturbations of these dynamos. 

We have refined the steady states from the branch with the resolution of 64^ x 32 Fourier 
harmonics. When solving the auxiliary problems and computing neutral eigenmodes of £* (the 
adjoint to the operator of linearisation), solutions in the form of truncated Fourier series have been 
sought. Pseudospectral methods [l^ have been used (resolution: 96^ x 48 Fourier harmonics before 
dealiasing) for computing various products encountered in the operator £. The code [l^ realising 
the biconjugate gradients stabilised method BiCGstab(f') j^-^| for £ = 2 has been applied to 
solve the linear system of equations for the Fourier coefficients. Table [T] shows the values of the 
coefficients in the system of amplitude equations for Ta = 675. 
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3.1. Linear stability of constant amplitudes 

The behaviour of solutions depends crucially on whether the evolutionary equations fl60l) - fl6T|) 
are parabolic. They constitute a parabolic (sub)system, when real parts of both eigenvalues of the 
matrix 

(dll Idi2 
/?21 /?22 


( 66 ) 


are positive. Otherwise, the second-order partial differential operator in the r.h.s. of fl 6 Up - (I^T]) 
would cause a superexponential growth of solutions in time (unless it is suppressed by nonlinear 
terms — but they are more likely to boost rather than inhibit this growth). The term 
is not involved in this analysis, because fl62|) is a non-evolutionary equation. The eigenvalues di 
shown in Table |2] reveal a high variability of the parameter values controlling the perturbation 
despite the interval, in Ta, of existence of the branch is short. At the left end of this interval for 
di > 0 but ^2 < 0, indicating a rapid (in the slow time) instability. In the remaining part of 
the interval, both dj > 0 and hence the second-order operator acts as the stabilising diffusion. 

Clearly, = 0 and any constant in time and space ^2 and ^3 (or, equivalently, C = 0 and 


Table 1: Coefficients in the amplitude equations for Ta = 675. 


Eq. 

Coeff. 

Values 

(EUD 

(dik 

0.884249123107103, 0.973164220817207, 0.002542972607168343 


7i 

0.588162048986537 

dSP) 

(d2k 

0.366161979311609, 4.15859432509866 


l2ij 

0.806971618405987, -4.23666924855939, 0.158691136470146, 
-4.82055714849942, -42.2469280336739, 0.135854143000776, 
0.209741278940542, 0.135854143000747, 0.003874598372548109, 


a 

0.192034252840662 

( 162 ]) 

iMj 

-7.08024457707350, 8.63396618933175, -0.500706409449314, 
8.63396618933175, 85.1114936402045, 0.385556218828258, 
-0.500706409449314, 0.385556218828260, -0.03382831817117658 

(ESI) 


1, -0.002659134450504667 


la 

0.596884895067222, -0.04846657396378542 


Table 2: Eigenvalues of matrix (IMl) . di, for two unit vectors q (such that are eigenvectors of matrix U). 


Ta 

q 

di 

d 2 

673.7 

qi = (0.998515, 0.054469) 
q2 = (0.822115, 0.569320) 

0.763328 

0.264885 

4.298197 

-261.047947 

675 

qi = (0.998613, 0.052645) 
q2 = (0.825155, 0.564905) 

0.778817 

0.270588 

4.264026 

-361.244316 

680 

qi = (0.998954, 0.045718) 
q2 = (0.836586, 0.547835) 

0.889724 

0.289251 

4.175111 

1326.420209 

690 

qi = (0.999477, 0.032314) 
q 2 = (0.858485, 0.512837) 

1.588269 

0.297058 

4.598703 

201.124913 

700 

qi = (0.999808, 0.019574) 
q 2 = (0.879679, 0.475566) 

2.518614 

0.174285 

16.785099 

212.326992 

703 

qi = (0.999873, 0.015897) 
q2 = (0.886016, 0.463653) 

2.597766 

0.057593 

65.565570 

281.662699 
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Figure 1: The region of linear stability (black) of constant solutions to the system of amplitude equations in the 
rectangle |co 3 | < 2, |co 4 | < 8, C = 0. 


any constant Cqs and C 04 ) are solntions to the system of amplitnde eqnations dSHD^dS^D, dSS])- In 
other words, any linear combination of the nentral modes Sfc _|_2 fl^TD with constant coefficients is a 
solntion to the stability problem in the broader class of weakly nonlinear large-scale pertnrbations. 
The condition of parabolicity of the linear operator in fl60|) - (l6T|) is the condition of stability of the 
solntion = 0 (i.e., (7 = Cqs = C04 = 0). A constant solntion (0,^2, ^3) is linearly stable if Rep < 0 
for both roots g of the qnadratic eqnation 


det 


/3ii + ^? 
1/^21 + r2i 
Fsi 


A 2 

1/^22 -b r22 + ip 
r32 


/^13 

r23 

-1 


0 


(see Fig. [T]). Here we have denoted 


3 

^2k = 

i =2 


^(73ifc + 73fci)^i I ( '^{'y3j3 + l33j)^j I - 


r3fc — 


0=2 


0=2 


1 + ( X](73i3 + 733i)0 
0=2 


3 

+ 'l3kj)^j 

i=2 


3.2. Two types of behaviour of amplitudes 

We report now results of numerical investigation of the system of amplitude equations 
([62]), fl65l) . Solutions, 27r-periodic in Y, have been sought as truncated Fourier series in the spatial 
variable. Pseudospectral methods [l^ have been applied. The presence of the cubic term in fl 6 Tl) 
necessitates performing the dealiasing by discarding a half of the highest-wave-number harmonics, 
rather than 1/3 of them as stipulated by the standard 3/2 rule 1^ (when evaluating quadratic 
terms). The second-order implicit midpoint method has been used to integrate the system in time. 
In terms of the Fourier coefficients of the variables ^s{y,T) denoted by ^(T) = (,^i(T ),..., 
the system takes the form 


dT 


= 1 < m < 2 M -b 1 , 


0 = fm{$.), 2M + 2<m<3M 


(67) 

( 68 ) 


(here evolutionary equations fl67|l stem from fl6UD . fIRT]) and fl6^ . and (1651) from (16^ ). In this 
notation, the midpoint method propagates solution at time tn a step forward to at time 

tn+i = tn + dt according to the equations 


tln+l) _ ^n) _ Uf f 
Sm Sm Jn 


^(n+l) _|_ ^(n) 


= 0, 1 < m < 2M -b 1, 


= 0, 2M + 2 < m < 3M. 


(69) 

(70) 
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To solve fl^ - dTO]) . we have used the routine newt [2l|] that realises a quasi-Newton method (which 
is a combination of the Newton method with a minimisation technique). The initial conditions are 
required to satisfy flTOjl . Mostly, a time step dt = 10“^ was used. A number of runs have been 
also performed using the standard fourth-order Runge-Kutta method and the time step dt = 10“^, 
with fITOl) being solved at each substep. Comparison of the results obtained by the two methods 
has shown a satisfactory quantitative agreement. 

Solutions to amplitude equations and fIBSl) apparently do not obey any conservation 

law and thus, in principle, can saturate at constant values or grow to infinity. In fact, we have not 
found any other patterns of behaviour. 

Let us now discuss in detail a run for Ta = 675 and q = (see Table [2]), which we will refer to as 
“the main sample run”. In a preparation of a blow-up, maxima of all the three amplitudes increase in 
time indefinitely; a graph of an amplitude near its maximum turns into a spike, travelling along the 
Y -axis at an increasing speed and hnally growing to very large values at a hnite time (see a typical 
behaviour in Fig. [2] and a video stored at arXiv.org as an ancillary hie Ta675.mp4). Simulated 
solutions were truncated to 256 Fourier harmonics in the spatial variable (i.e., they involved wave 
numbers up to 127, and computation of nonlinear terms involved 512 harmonics before dealiasing). 
The noise in the graphs of absolute values of the Fourier coefficients of the amplitudes (panels (j)~(l) 
in Fig. 12]) at T = 31 for wave numbers exceeding 95 indicates, that such a resolution is excessive 
when computing with the double precision up to this time. In the course of evolution, the energy 
gets eventually almost evenly distributed between all Fourier harmonics (panels (j)-(l) in Fig. |2|); 
at the hnal stage, amplitude C resembles a periodically-replicated ^-function moving at a constant 
speed. (Of course, at this stage the numerical solution is underresolved and cannot be fully trusted 
any more.) A common feature of nonlinear partial differential equations of the hrst order such as 
the inviscid Burgers equation is development of a jump due to intersection of characteristics; in the 
present case the singularity can be rather interpreted as a derivative of a function, experiencing a 
jump. However, addition of a small diffusivity into the Burgers equation smooths out the jumps, 
which is not the case in the present system of amplitude equations, where diffusion is present and 
is not inhnitesimal. 

Blow-ups are well-known to be linked with self-similar solutions (see, e.g., 0 , 0 )- It is easy to 
check that equations (I60|) - (l62|) and ( 16 ^ admit self-similar solutions of the form 

(71) 


(C,c„3,c„4) = (T. - r)-‘/"<*(y(r. - 


(see 0 ). To examine, whether the solution has this asymptotics near the singularity, we plot in 
Fig. El graphs of the quantities 


Mf,raUT) = l/max/(F,T), M^,^i„(T) = l/mm/(F,T), (72) 

as well Af|j^g^^(T) and M|j^j^(T); here the function / is any amplitude C, Cqs and C 04 . (The maxima 
and their positions are determined applying the quadratic interpolation.) The graphs reveal that 
the factor (T* — may indeed correctly describe the growth of the amplitudes at intermediate 

times preceding the blow-up (31 < T <33 for the maxima of the solution under consideration, as 
seen in panels (g)-(i) of Fig. |3|). (This is not seen as clearly in panels (j)-(l) of Fig. [3]for the minima, 
because miny co 4 (R, T) vanishes at T 33.66, forcing us to plot Mjmin(T) = 1/ miny/(F, T) for 
a shorter time interval.) However, further on, at times immediately preceding the blow-up, the 
maxima and minima behave as (T* — T)~^ (as follows from graphs in panels (a)-(f) of Fig. |3|). 
The graphs in panels (a)-(c) for the maxima of amplitudes and in panels (d)-(f) for their minima 
coherently indicate that the blow-up occurs at T* 34.17. 

To see how well the self-similar solution (ITT]) approximates near the singularity the spatial de¬ 
pendence of the amplitudes, we have explored the behaviour of Pf{T), the position of the maximum 
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Figure 3: The quantities ^ (vertical axis) M/,max(T) (panels (a)-(c)), Mf^^iniT) ((d)-(f)), ((g)- 

(i)) and ((j)-(l)) for C ((a), (d), (g), (j)), Cqs ((b), (e), (h), (k)) and Co 4 ((c), (f), (i), (1)) during the 

preparation of the blow-up of a solution to the system of amplitude equations. Horizontal axis: the slow time T. 
Dashed lines on panels (a)-(f): lines through the solid points on the respective graphs, approximating the asymptotes 
that characterise the final stage of the process of preparation. 
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Figure 4: Positions of the maxima (vertical axis), Qf(T), (panels (a)-(c)) and Q‘f{T) ((d)-(f)) for C ((a), (d)), cqs 
(( b), (e)) and co 4 ((c), (f)) during the preparation of the blow-up of a solution to the system of amplitude equations. 
Horizontal axis: the slow time T. Dashed lines: same quantities computed using a double spatial resolution and a 
10 times smaller time step. 

of f{Y,T) at time T, where as function / we again take the amplitudes C, Cqs and C 04 . For the 
self-similar ansatz (ITT]) we would obtain 

Pf{T) = P^{T,-Ty/\ (73) 

where P^ is the argument at which the prohle $ in (ITTll takes the maximum value; we therefore 
expect Pf{T) to be a linear function of time. Graphs of positions of the maxima, Pf{T), as well 
as Pf(T) as functions of the slow time T are plotted in Fig. 01 and they do not contradict the 
hypothesis, that asymptotically d73|) holds in the time interval 33 < T < 34. However, closer to 
the time of singularity T*, the graphs of Pj{T) (panels (d)-(f) of Fig. 0]) bend down. To check that 
this bending is genuine and not caused by the lack of resolution, we have recomputed the solution 
with a twice higher spatial resolution and a 10 times smaller time step; the graphs obtained for 
the rehned positions are shown in Fig. 0]by dashed lines, which visually coincide with the original 
graphs. Note, that the graphs of Pf{T) (panels (a)-(c) of Fig. 0]) feature a similar behaviour, which 
is close to a linear one at a larger time interval; however, Pf{T) bends down near T* more visibly. 
Comparison of panels (a)-(c) versus panels (d)-(f) shows, that it is difficult to make a reliable 
conclusion about the exponent in the argument of $ in (1731) based on this data. 

We have investigated how the occurrence and the time of blow-up depend on the mean am¬ 
plitudes of the initial conditions (see Fig. |5]). The mean amplitudes ((cqs)) and ((C 04 )) of the two 
neutral zero-mean small-scale stability modes, that exist due to translation invariance in horizontal 
directions of equations of convective dynamo, appear to be important control parameters: if in the 
course of temporal evolution these means quit the region of linear stability of constant solutions 
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Figure 5: Isolines of blow-up times for varying initial ((coa)) and ((C 04 )). The region, for which solutions tend to 
constant values, is cross-hatched. When a point (((coa)), ( 004 }) approaches this region, the times of blow-up tend to 
infinity, and hence there are isolines for the values 10 and higher near the upper border of this region, and for values 
1 and higher near the lower border. These isolines are very close to the respective borders and are not shown. The 
region of linear stability of constant solutions to the system of amplitude equations is shown in gray (cf. Fig.[T]). 


(see the previous subsection), typically, the solution eventually experiences a blow-up at a hnite 
time. 

Ideally, in such an investigation we would like to compare solutions to the system and 

( 16 ^ obtained for initial conditions where these means are varied and the fluctuating parts are kept 
unaltered. However, this is impossible, because the initial data must satisfy the non-evolutionary 
equation ([62j). Consequently, Fig. |5]is plotted for solutions with the initial conditions constructed 
as follows: At T = 0, the fluctuating parts of the amplitudes Co 3 (H) and Co 4 (V) are the same in all 
the runs (they coincide, respectively, with the fluctuating parts of Cqs and C 04 in the discussed above 
main sample run at T = 30; the r.m.s. means of the fluctuating parts are, respectively, 0.22 and 
1.88). The zero-mean amplitude C(V) is then found from fl62|) for each pair ((cqs)), ((C04)) separately. 

We see in Fig.[5]that the region of mean amplitudes, for which solutions tend to constant values, 
is within the region of linear stability of constant solutions (see subsection 13.11) . We also observe 
that when a point (((cqs)), ((C04))) approaches the lower boundary of this region outside it, the times 
of blow-up tend to infinity much faster than when the upper boundary is approached, suggesting 
that the blow-up preparation proceeds differently, at least in details, in the two areas above and 
below the region of linear stability of constant solutions. Simulations support this conclusion: For 
the main sample run, ((cqs)) = 0.45, ((C 04 )) = 8.02 at T = 30; it corresponds to a point in the area 
above the region of stability of constant solutions in Fig. O An additional run for Ta = 675 and the 
initial conditions, that have been constructed as explained above for ((cqs)) = — 1.2 and ((C 04 )) = 6 , 
is represented in Fig. [5] by a point in the area below the region (a video is stored at arXiv.org as 
an ancillary hie Ta675a.mp4). In both runs the graphs of the moduli of Fourier coefficients of the 
amplitudes versus the wave number k remain near the blow-up predominantly linear. However, 
while in the main sample run their steepness monotonically decreases in time, in the additional 
run the steepness experiences a complex oscillatory behaviour. Also, while at the hnal stage one 
spike of a large height develops near the time of singularity in the main sample run, high-amplitude 
spatial oscillations burst out in space in the additional run. 

S.S. Dependence on the Taylor number 

We summarise now the numerical results obtained for several other steady states from the 
branch Although the coefficients of the amplitude equations signihcantly vary within the 

branch (most of them by more than an order of magnitude), we do not encounter any other types 
of the temporal behaviour of amplitudes, rather than those found for Ta = 675: the development 
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Figure 6: Lebesgue norms (vertical axis) of amplitudes cqs (upper panel), co4 (middle panel) and C (lower panel) 
as functions of the slow time T (horizontal axis; the scale is logarithmic to have a better display of details of the 
graphs for Ta > 690, for which the blow-up times are small) for q = (see Table [2]). The plots are labelled by the 
respective Taylor number values. 

of a singularity at a finite time, or the decay of coefficients to constant values, i.e., degeneration 
of a large-scale perturbation to a neutral zero-mean short-scale mode existing due to the spatial 
invariance of the convective system. 

The genericity of these two types of behaviour is illustrated by Fig. |6], showing the temporal 
evolution of the Lebesgue norms of amplitudes (defined as for an amplitude c{Y,T)) 

for several values of the Taylor number. The graphs in Fig. |6]are plotted for the initial conditions 
constructed essentially by the same procedure as in the previous subsection: At T = 0, amplitudes 
co 3 (F") and co 4 (F") are the same in the runs for all considered Ta, and coincide with those in the 
main sample run at T = 30. The zero-mean amplitude C{Y) is found from (l62|) for each Taylor 
number separately. 

Figure E] shows that the blow-up times depend on the Taylor number non-monotonically. Al¬ 
though for Ta > 690 the employed initial profiles of Cqs and Co 4 quickly lead to a blow-up, for the 
nearby value Ta = 680 the singularity does not develop and the perturbation tends to a short-scale 
zero-mean neutral mode. Noteworthily, prior to the hnal stage amplitude C remains close to unity. 
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and therefore the cubic nonlinearity of the system of amplitude equations, which is associated with 
= C (see dSID), apparently does not play an important role in the development of the singularity. 

Further qualitative information is revealed by the videos of the hnal stage of the preparation of 
the blow-up in two runs for Ta = 680 and 703 (for initial conditions constructed by an alternative 
procedure), and in the additional run for Ta = 675 (in all the three runs q = q^), which are stored 
at arXiv.org as ancillary hies Ta680.mp4, Ta703.mp4 and Ta675a.mp4, respectively. For Ta = 703, 
like in the runs for Ta = 675, the graphs of the moduli of Fourier coefficients of the amplitudes as 
functions of the wave number k remain predominantly linear, although the hne structure of these 
graphs is much richer than in the case of the main sample run — all amplitudes have now several 
extrema, whose magnitudes grow in time. In the main sample run for Ta = 675, the steepness of 
such graphs just decreases in time; by contrast, for Ta = 703 initially the steepness increases, so that 
the round-off level ~ 10“^® is attained at A; ~ 60 (in this simulation, solutions have been truncated 
to 512 Fourier harmonics in Y after dealiasing), and afterwards the slope starts to decrease; in the 
additional run for Ta = 675 the slope experiences an oscillatory behaviour. Eventually the slope 
of the graphs decreases until near the time of singularity the graphs become close to horizontal 
ones. The growth of the extrema (in absolute values) of the amplitudes intensihes in time, but 
none of the extrema gives rise to a ^-function, as the single maximum of C does in the main 
sample run for Ta = 675. For Ta = 703, just before the emergence of the singularity, a secondary 
fast moving perturbation of amplitudes develops, which modihes their so far preserved forms; the 
overall linear shape of the graphs of absolute values of the Fourier coefficients becomes deformed. 
In the additional run for Ta = 675, high-frequency high-amplitude spatial oscillations are observed. 
However, the spatial resolution becoming insufficient, the occurrence of these hnal events cannot 
be fully trusted. 

The evolution of perturbations near the singularity for Ta = 680 signihcantly differs from those 
for Ta = 675 and 703. The singularity sets in much faster. Two and then three global maxima 
develop in graphs of the moduli of Fourier coefficients of the amplitudes versus the wave number, 
which are initially predominantly linear. Each amplitude has initially one maximum; the maxima do 
not yield (5-functions, but rather burst into packets of amplitude-modulated blobs of high-frequency 
oscillations. 

4. Concluding remarks 

We have performed a two-scale analysis of large-scale perturbations of a periodic array of steady 
thermal convective hydromagnetic states constituting the branch j^. These short-scale dy¬ 
namos are symmetric about the vertical axis and stable to short-scale perturbations. We have 
shown that, in a collective action, they are capable of generating a large-scale magnetic held, whose 
geometry is determined by the properties of the short-scale operator of linearisation around the 
respective steady state. The operation of this large-scale dynamo is based on the joint action of 
the combined eddy diffusivity and eddy advection, but it does not involve the (a-effect. To the best 
of our knowledge, this mechanism is new. 

Several remarks are in order. 

1°. We have found that only two small-scale neutral stability modes of a steady convective 
dynamo with such a group of symmetries have non-vanishing means, each involving a combination 
of a mean horizontal velocity and magnetic held. The solenoidality in slow variables of the two mean 
helds implies that only a specihc linear combination of these two modes is involved in the asymptotic 
expansion of a large-scale perturbation, depending on a slow spatial variable and evolving in the 
slow time. For a given state W from the considered branch, two such linear combinations can be 
determined. Two other neutral small-scale stability modes, d'W/dxi, existing due to translation 
invariance in horizontal directions of equations of convective dynamo, are zero-mean, i.e., their 
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amplitudes are not associated with any mean fields. The two amplitudes are essential for the large- 
scale dynamics of the perturbation: if they are both set to zero, then, by virtue of equations ( 1 ^ . the 
amplitude of the remaining small-scale neutral mode, involving non-zero mean-field components, 
also vanishes. 

This clearly shows that omitting amplitudes of neutral zero-mean small-scale stability modes 
(when they exist) in the description of the dynamics of multiscale perturbations, which is the com¬ 
mon practice in the theory of mean-field electrodynamics, can render the description inadequate. 
Such modes are not rare; they are featured by translation-invariant and/or time-invariant MHD 
systems, and can be obtained by differentiating the states, whose stability is analysed, in the spatial 
variable along the invariant direction and/or time. 

2°. Requirements for the averaging procedures that can be applied also differ significantly be¬ 
tween MFE and MST. In MFE, any averaging is deemed acceptable when analysing MHD turbu¬ 
lence, provided it satisfies the Reynolds rules. Accordingly, averaging over a horizontal periodicity 
cell (or, equivalently, over entire horizontal planes) is often considered in literature. Perhaps, such 
averaging could be justified in the framework of statistical approach to turbulence (e.g., for en¬ 
sembles of turbulent eddies). By contrast, our derivation of amplitude equations, following the 
MST approach and using the PDE homogenisation techniques, shows that the only appropriate 
averaging when considering two-scale perturbations of laminar MHD regimes is over a periodicity 
domain of the perturbed state; for time-periodic perturbed states, this must be supplemented by 
averaging over the temporal period. 

Thus, a mathematically rigorous asymptotic procedure leaves no freedom for choosing averaging 
to our taste; using spatial (or spatio-temporal) averaging is inevitable — this follows from the 
requirement of solvability of equations in fast variables. Consequently, the MFE approach to such 
problems (e.g., kinematic dynamo problems concerning generation of large-scale magnetic field by 
a periodic array of small-scale flow cells) may yield quantitatively incorrect results, when other 
types of averaging are used, even when they satisfy the Reynolds rules and are allowed in MFE. 
Furthermore, no conclusions of potential interest for astrophysics can be drawn from the results 
of 0 , 1 ^ on the significant (i.e., responsible for the evolution of large-scale magnetic field) a-effect, 
because averaging over half a periodicity cell performed in these works does not satisfy the Reynolds 
rules (the significant a-effect ibid, is zero). 

3°. The system of amplitude equations that we have derived and studied numerically involves 
terms describing eddy diffusivity and eddy advection, but not the a-effect. This calls for a com¬ 
parison of the action of these mechanisms of magnetic field generation important in astrophysics. 

The spectrum of the a-effect operator H i—)■ V x 21H, acting, for instance, on space-periodic 
fields H, is symmetric about the imaginary axis (see also j^, section 3.3.2). Under the action of 


this operator, a generic perturbation experiences a superexponential growth in slow time, which in 
this case is T = et (while in the absence of the a-effect large-scale amplitudes evolve in 0(e) slower 
time T = e‘^t, see section [53]). In the MST framework, when 21 7 ^ 0, the a-effect operator typically 
controls the evolution of mean fields of perturbation solely (see [33|, sections 3.3.2, 4.2.3, 6.3.2, 
7.3.2, 8.3.2). The superexponential growth might be restrained by diffusion, since it is governed 
by a higher-order differential operator. However, the orders of the a-effect and eddy diffusivity 
differential operators being different, the operators arise in solvability conditions at different orders 
of the scale ratio and do not coexist. In principle, emergence of the operators of a-effect and eddy 
diffusivity together is not ruled out in the MFE theory, where the so-called /9-effect is sometimes 
considered, in whose presence the mean e.m.f. becomes 


(V X H) = 5 ^ 






Up./''’ 

i,j,k 


'dX, 




However, this equality is only partially justified in the MFE theory by considering the truncated 
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Taylor expansion of the mean field in the integral operator relating the flnctnating part of the 
magnetic held to the mean part, and assnming that the respective kernels decay fast (see 22|). 


Whether this jnstihcation can be extended to a formal asymptotic argnment is an interesting open 
qnestion. 

Thns, to the best of onr knowledge there are no mathematically solid evidences that in mnltiscale 
problems, where a rigorons derivation of amplitnde eqnations is possible, the action of the a-effect 
operator can be restrained by eddy diffnsivity or the mechanisms of nonlinear satnration. Hence, it 
can be described as essentially “self-destrncting”: In a system with the a-effect, the pertnrbation 
grows snperexponentially nntil it loses the asymptotic smallness and thereby signihcantly affects 
the pertnrbed state. In other words, the a-effect can generate magnetic held for a relatively short 
time, bnt it destabilises the existing small-scale strnctnre and modihes the entire MHD system 
into a completely new regime. For this new regime, the analysis of eddy ehects must be done 
anew (provided the scale separation is preserved and the multiscale analysis can be implemented; 
actually, it is also not guaranteed that the new regime will be capable of generating the held). 
This is incompatible with the astrophysical reality — large-scale cosmic magnetic helds are known 
to exist for very long times. Therefore, it is natural to expect that a sequence of such events of 
self-destruction terminates at an MHD state, where the magnitude of the a-ehect is sufficiently 
small not to cause further disturbances. 

Within this line of reasoning, the feasibility of the a-ehect paradoxically, may owe to the phe¬ 
nomenon, originally perceived as signihcantly reducing the importance of the a-ehect, namely, the 
a-quenching, i.e., a drastic decrease of the a-ehect on increasing the magnetic Reynolds number, 
discovered by S.I. Vainshtein [^. In the kinematic dynamo problem for a how with an external 
scaling, the a-ehect operator is present in the equation for mean helds together with molecular 
dihusion (see also 


33| . ch. 10, 11). In astrophysics, molecular dihusion is weak. However, 


if the a-quenching inhibits the action of the a-ehect to the levels, where it is essentially ohset 
by molecular magnetic dihusivity, then, at least within this model (also mathematically rigorous), 
the a-ehect operator may contribute to stationary magnetic held generation without causing the 
restructuring of the MHD system hosting it. 

As our results show (see section [3]), in the absence of the a-ehect the overall character of the 
large-scale evolution is not that diherent. If negative combined eddy dihusivity is present, then a 
generic perturbation grows snperexponentially. As mentioned above, the characteristic slow time of 
the large-scale evolution is then slower than for the a-ehect, but the ultimate result of the evolution 
is the same: destruction of the underlying MHD regime and the onset of a new one (which can lack 
scale separation). We have no evidence that taking into account nonlinear terms in equations for 
weakly nonlinear perturbations can halt such processes. Here we have studied numerically the case, 
where the action of the combined eddy dihusivity is stabilising (as that of molecular dihusivity). 
In this case the large-scale evolution signihcantly depends on how large is the initial perturbation. 
When it is below a certain threshold, the amplitudes tend to constants and the perturbation evolves 
to a linear combination of the two zero-mean neutral modes; otherwise, it blows up in apparently 
a hnite time. Therefore, like in the presence of the a-ehect, the large-scale evolution controlled by 
eddy dihusivity and eddy advection can perturb the original small-scale dynamics till its structure 
is signihcantly modihed, but the growth rate of these processes is much, 0(e), slower. 

An open question is whether the system of amplitude equations and has a solution 

which is a travelling wave. We have failed to hnd this regime for the steady dynamo at Ta = 675, 
but this can be due to the scarceness of our ehorts; also, it can, in principle, exist for perturbations 
of other steady states from the branch that we have considered, or in other branches with the 
symmetries compatible with our multiscale analysis (such as the time-periodic dynamo Pj^^). If 
such a solution exists and is stable, it would be an example of a large-scale perturbation that does 
not act self-destructively as just discussed; if it is unstable in slow variables, it is nevertheless of 
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interest as an invariant object in the phase space influencing the dynamics of amplitudes. 

4°. We have analysed numerically the system of amplitude equations derived here for convective 
hydromagnetic regimes with the symmetry group of states in the branch The behaviour of 
its solutions is complex, because this is a mixed system, involving both evolutionary and non¬ 
evolutionary nonlinear equations. Many questions remain open, for instance: Does this system 
allow patterns of behaviour of perturbation, different from those that we have identified? Does the 
blow-up occur for other boundary conditions? Does the solution for Ta = 675 “converge” at the 
blow-up time to a distribution involving the ^-function, as computations suggest? If this ^-function 
is interpreted as a derivative of certain quantities experiencing a jump, can these quantities be 
identified in physical terms? How can a solution be continued in some weak sense beyond the 
blow-up time T*? 

5°. The system of amplitude equations that we have derived is very restrictive: it can only be 
used to describe weakly nonlinear perturbations, where amplitudes are initially constant on lines, 
perpendicular to certain well-defined directions q on the plane of slow variables. This system, 
therefore, does not describe general two-scale perturbations of a convective MHD state, that are of 
the form flT^ - ffTTl) and obey equations fl8|)- ffTT|) . An interesting open question is to construct an 
asymptotic formalism for perturbations of a more general form. 

6 °. Finally, it is of interest to perform similar investigations for other branches of short-scale 
convective hydromagnetic regimes: (i) for time-periodic MHD states constituting branch that 
emerge in the interval 704 < Ta < 705 in a Hopf bifurcation from considered here (large- 
scale perturbations of time-periodic states involve an additional amplitude); (ii) for steady states 
and time-periodic regimes constituting the branches, found in j^, that possess other groups of 
symmetries, for which 77 = 0. 
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